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Purpose: Eosinophils function as an effector cell in the development of asthma and allergic disease. Eotaxins are cytokines that promote pulmo- 
nary eosinophilia via the receptor CCR3. Single-nucleotide polymorphisms (SNPs) in CCR3 and eotaxin genes are associated with asthma. In this 
study, genetic interactions among SNPs of several eotaxin genes and CCR3 were assessed and their relationship with blood eosinophilia in asthma 
was examined. Methods: A total of 533 asthmatics were enrolled in this study. Asthmatics with eosinophilia (> 0.5 x10 a /l_) were compared with 
those without eosinophilia (<0.5x10 9 /l_). Chi-square tests were used to compare SNP frequencies. Two different models were used to evaluate 
gene-gene interactions: logistic regression and generalized multifactor dimensionality reduction (GMDR). Results: E0T2+304C> A{2°L>\) was 
significantly associated with 3 of the 4 CCR3 SNPs among asthmatics with eosinophilia (P=0.037-0.009). E0T2+304O A (22l>\) and the CCR3 SNPs 
were also significantly associated with blood eosinophilia in an interaction model constructed by logistic regression (P=0.0087). GMDR analysis 
showed that the combination of E0T2+304OA (29L>I) and CCR3- 174C> 7" was the best model (accuracy=0.536, P=u.005, CVC 9/1 0). Conclu- 
sions: The epistatic influence of CCR3 on eotaxin gene variants indicates that these variants may be candidate markers for eosinophilia in asthma. 
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INTRODUCTION 

Eosinophils play an important role in the development of asth- 
ma and allergic disease as an effector cell. The extent of eosino- 
philic inflammation is a determinant of the severity of asthma 
symptoms, 1 and it usually correlates with airflow Limitations. 2 
Relationships between eosinophilic infiltration and Th2 cyto- 
kines have been identified within the target tissues of allergic 
disorders. 3 The processes of eosinophilic infiltration and periph- 
eral eosinophilia are dependent on eosinophil-specific cyto- 
kines and chemokines {e.g., IL-5; the eotaxin family; RANTES; 
and MCP-2, -3, and -4), which provoke an eosinophilic response 
in the peripheral blood and airways via CCR3. 4 ' 5 CCR3 mRNA 
and protein levels are elevated in the bronchial mucosa of asth- 
matics, and this elevation is associated with airway hyperre- 
sponsiveness. 6 The participation of CCR3 in airway eosinophil- 
ic infiltration has been demonstrated in a study of CCR3-defi- 
cient mice, which showed that most eosinophils are arrested in 
the subendothelial space. 7 Eosinophils are derived from CD34(+) 
hematopoietic progenitor cells. 8 Eotaxin- 1 and -2 induce the 
migration of bone marrow and blood CD34(+) CCR3(+) cells in 



vitro? Therefore, the CCR3-eotaxin pathway is important in the 
regulation of allergen-induced hematopoiesis and in the accu- 
mulation and mobilization of eosinophil lineage-committed 
progenitor cells in the lung. 9 The human CCR3 gene (MIM 
#601268) is located on chromosome 3p21.3. 10 The single-nucle- 
otide polymorphisms (SNPs) CCR3-22557G>A and -1740T 
are associated with the number of eosinophils in asthmatic pa- 
tients. It has previously been reported that CCR3 expression is 
higher on the eosinophils of asthmatic patients lacking a haplo- 
type composed of rare alleles of 4 CCR3 SNPs. 11 

Eotaxin- 1 is important for eosinophilic inflammation in early 
stages of the asthmatic response, while eotaxin-3 may account 
for eosinophil recruitment to the airways in late stages of the 
asthmatic response. 12 Three eotaxin family members, eotaxin- 1, 
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-2, and -3, are selectively bound by CCR3. Eotaxin-1 and -2 show 
different pathophysiologic responses in relation to eosinophils. 
Disruption of the eotaxin-1 gene leads to a slight reduction in 
the eosinophil count in the blood and airways. 13 Eotaxin-2-de- 
ficient mice have normal baseline eosinophil levels in their tar- 
get tissues, and they do not develop airway eosinophilia in re- 
sponse to IL-13. 14 Two eotaxin polymorphisms, EOT2+1272A>G 
and E0T1+123G>A, are associated with asthma and high serum 
total IgE levels, respectively. 15 

Because asthma is a multifactorial disease, the genetic com- 
ponent may be derived from the combined effect of numerous 
genes. Individual genes may act independentiy or in combina- 
tion with other genes in the same biological pathway, resulting 
in variable effects. 16 The CCR3 gene and eotaxin gene family are 
possible contributors to the development of asthma or other al- 
lergic phenotypes via the receptor-ligand interaction. 

In this study, we evaluated gene-gene interactions between 
CCR3 and eotaxin in relation to eosinophilia in patients with 
asthma. 

MATERIALS AND METHODS 
Subjects 

We enrolled 533 Korean asthmatics in our study. The Institu- 
tional Review Board of [Soonchunhyang University Bucheon 
Hospital] approved the study. All of the patients had current 
symptoms, including wheezing, dyspnea, and cough, and met 
the criteria for asthma as defined by the American Thoracic So- 
ciety. 17 Each patient showed airway reversibility as document- 
ed by an inhalant bronchodilator-induced improvement in the 
forced expiratory volume in one second (FEV1) of > 15% 18 and/ 
or airway hyperresponsiveness as shown by a provocative con- 
centration of methacholine required to cause a 20% decrease in 
the FEV1 (PC20) of <8 mg/mL. 18 The asthma subjects were di- 
vided into 2 groups according to the presence of eosinophilia. 
Eosinophilia was defined as a blood eosinophil count exceed- 
ing 0.5 X 10 9 /L. 19 No subject had used systemic or inhaled ste- 
roids for 4 weeks before peripheral blood eosinophil counting. 
No subject had a history of medication or the ingestion of raw 
foods that could produce blood eosinophilia. ELISAs for 4 com- 
mon parasites {Paragonimus, Clonorchis, Cysticercus, and Spar- 
ganum) were negative in the enrolled patients. The clinical 
characteristics of the subjects are summarized in Table 1. 

Genotyping by single-base extension and electrophoresis 

To genotype polymorphic sites, amplifying and extension 
primers were designed for single-base extension. All primer ex- 
tension reactions were performed using a SNaPshot ddNTP 
Primer Extension Kit (Applied Biosystems, Foster City, CA) ac- 
cording to the manufacturer's instructions. Detailed methods 
are described in the online supplement. 



Table 1. Clinical characteristics of asthma subjects according to the presence 
of eosinophilia 





Asthma wvith 

AAOLI IIIIG Will 1 

eosinophilia 


Asthma \A/rthni it 

rAO LI 1 1 1 1 a V V 1 LI 1 U LI L 

eosinophilia 


Pvalue 


No. of subjects 


155 


378 




Median age (range) 


37 (8-75) 


43 (8-81 ) 


<0.0001 


Sex (males/females) 


92/63 


1 52/226 


<0.0001 


Smokers (%) 


20.4 


23.8 


0.42 


rvL /o, preo. 


00 Q _l_ 1 A 

oo.o 3: 1 A 


qq n-t-n q 
oo. U 3: u.y 


U.oo 


FEV1%,pred. 


81.1 ±1.8 


80.8 + 1.1 


0.86 


% change in FEV1 by bronchodi- 
lator 


11.5±1.4 


10.5 + 0.6 


0.54 


PC20 (mg/mL) 


1.7 + 0.2 


3.0 + 0.3 


0.0001 


Total IgE (lU/mL) 


560 + 64.1 


349+27.2 


0.003 


Blood eosinophil count ( x 1 0 3 /L) 


0.891+0.04 


0.223 + 0.01 


<0.0001 


Blood eosinophil count (%) 


10.9 + 0.4 


3.4+0.1 


<0.0001 



Eosinophilia was defined as a blood eosinophil count >0.5X10 9 /L. P values 
were calculated for comparisons between asthmatics with eosinophilia and 
without eosinophilia. Chi-square tests were used to compare categorical vari- 
ables; Student's f-test was used to compare continuous variables. Mean ± SEM. 



Analysis of polymorphisms, Hardy-Weinberg equilibrium, and 
linkage disequilibrium 

A total of 14 SNPs in eotaxin-1, -2, and -3 (6 SNPs in eotaxfn-1, 
5 in eotaxin-2, and 3 in eotaxin-3), and 4 SNPs in CCR3 were in- 
cluded in our analysis. Testing for Hardy-Weinberg equilibrium 
and the calculation of D' for the identification of linkage disequi- 
librium were performed using PHASE v2.0.2 20 and Arelquin 
v2.0. 21 Continuous variables, such as the blood eosinophil count 
and percentage, serum total IgE level, and PC20 methacholine 
level, were log-transformed to approximate normal distribu- 
tions. Differences in the log-transformed values among geno- 
types were examined using a generalized linear model type III 
SS. Statistical significance was defined at the standard 5% level. 

Dependency of allele frequency and logistic regression 
interaction modeling 

The dependency between each pair of CCR3 and eotaxin 
SNPs in asthmatics with and without eosinophilia was tested by 
Chi-square analysis. 

An epistatic model was designed based on the logistic regres- 
sion model. 22 "CCR3e£f," a new numerical variable as an effect of 
CCR3 SNPs on blood eosinophilia, was derived according to the 
extent of contribution to eosinophilia in asthma using logistic 
regression. Odds ratios (ORs) for eosinophilia for each geno- 
type at each locus were calculated in asthmatics with eosino- 
philia. Each genotype was coded according to the ORs as fol- 
lows: homozygote with the highest OR, 2; heterozygote, 1; and 
the other homozygote, 0. 

In CCR3 ef f=ai {-22557G> A)+a 2 (-520T>G)+a 3 {-174C>T)+a A 

(+51T>C), 
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a denotes the coefficient fitted by logistic regression (ai =0.3601, 
a 2 =0.3525, a 3 =0.4204, and a 4 =0.1867). Epistasis was subse- 
quently tested independently for each SNP in eotaxin by fitting 
the following logistic regression model: 

P (Eosinophilia)=exp (X)/[l+exp (X)] 

where X=/U/?iCCR3 e ff-^2SNP E ot™+A (CCR3effXSNPEot™n), fio 
is the intercept, ji represents the parameter estimates, and SNP 
Eotaxin IS each eotaxin genotype (0, 1, or 2). 

The Bonferroni correction for multiple comparisons was used 
for each eotaxin gene (the global significance level was adjusted 
to P<0.0083 for 6 eotaxin-1 SNPs, P<0.01 for 5 eotaxin-2 SNPs, 
and P<0.017 for 3 eotaxin-3 SNPs). Statistical significance was 
defined at the standard 5% level. 

Generalized multifactor dimensionality reduction (GMDR) 

To validate the results of our logistic model, we performed 
GMDR analysis. First, the data were randomly split into 10 equal 
parts for cross-validation. Nine of these were used as a training 
set; the remaining set was used for independent testing. Cross- 
validation is a measure of the number of times a particular set 
of loci is identified in each possible 9/10 of the subjects. Second, 
a set of n genetic factors was selected. Third, all possible multi- 
factor classes or cells were represented in rc-dimensional space 
and a cumulative score was calculated within each cell. 23 Fourth, 
each multifactor cell in n- dimensional space was labeled as 
"high-risk" if the average score met or exceeded the threshold 
(0) or as "low- risk" if the threshold was not exceeded. This pro- 
cess was repeated for each possible cross-validation interval. 
Fifth, all potential combinations of n factors were evaluated se- 
quentially for their ability to classify cases and controls in the 
training data, and the best rc-factor model that yielded the min- 



imum misclassification error was chosen. Sixth, the indepen- 
dent testing set was used to estimate the prediction error of the 
best model selected in the fifth step. Finally, among this set of 
the best models, we pick the model with rninimum prediction 
error and/or maximum cross-validation consistency 23,24 A co- 
variate analysis was performed as described above 23 with age 
(continuous variable) and sex (discrete variable) as covariates. 

RESULTS 

The genotype distributions of the four SNPs in CCR3 and 14 
SNPs in the eotaxin genes were in the Hardy- Weinberg equilib- 
rium (data not shown). No pair of SNPs showed strong linkage 
disequilibrium ( > 0.8) between the eotaxin gene SNPs and CCR3 
SNPS (see Supplementary Table S2). 

Epistasis was evaluated to access whether the eotaxin SNPs 
modified the effect of the CCR3 SNPs on blood eosinophilia 
susceptibility. We examined whether dependence existed be- 
tween the CCR3 SNPs and eotaxin gene SNPs before the epista- 
sis tests. E0T2+304OA was dependent on 3 intronic SNPs 
(CCR3-22557G>A, CCR3-520T> G, and CCR3-1 74C> T) in CCR3 
among asthmatics with eosinophilia by Chi-square analysis 
(Table 2). 

The epistatic model constructed by logistic regression that 
tested the effect of the combined CCR3 SNPs on blood eosino- 
philia differed depending on the eotaxin SNP: E0T2+304OA 
(29L>I) with the CCR3 SNPs was significandy associated with 
blood eosinophilia (P= 0.0087; Table 3). Thus, genetic suscepti- 
bility to blood eosinophilia in asthmatics is conferred by both 
CCR3 and eotaxin. 

According to our GMDR results, the combination of CCR3- 



Table 2. Tests for dependence between the CCR3 gene SNPs and eotaxin gene SNPs 

With eosinophilia Without eosinophilia 

Eotaxin CCR3 — 







-22557G>A 


-520T>G 


-1740T 


+51T>C 


-22557G>A 


-520T>G 


-1740 T 


+51T>C 


Eotaxin-1 


-1329G>A 


0.13 


0.37 


0.91 


0.97 


0.41 


0.24 


0.043 


0.45 




-329A>G 


0.43 


0.71 


0.8 


0.84 


0.54 


0.84 


0.95 


0.49 




+123G>A 


0.61 


0.93 


0.77 


0.17 


0.12 


0.042 


0.47 


0.23 




+371 T> A 


0.17 


0.52 


0.85 


0.5 


0.40 


0.84 


0.51 


0.8 




+1B54G>A 


0.3 


0.29 


0.77 


0.83 


0.83 


0.43 


0.43 


0.69 




+1671T>C 


0.99 


0.037 


0.26 


0.62 


0.66 


0.71 


0.87 


0.55 


Eotaxin-2 


+9BG>A 


0.83 


0.75 


0.67 


0.92 


0.33 


0.12 


0.41 


0.61 




+304OA 


0.037 


0.016 


0.009 


0.53 


0.65 


0.74 


0.41 


0.54 




+4470 T 


0.85 


0.85 


0.69 


0.37 


0.84 


0.92 


0.58 


0.14 




+1272 A>G 


0.64 


0.16 


0.09 


0.45 


0.77 


0.59 


0.39 


0.86 




+1923A>C 


0.91 


0.57 


0.83 


0.47 


0.87 


0.36 


0.4 


0.77 


Eotaxin-3 


+77C> T 


0.24 


0.94 


0.94 


0.95 


0.83 


0.39 


0.39 


0.15 




+71B4A>G 


0.33 


0.2 


0.1 


0.52 


0.42 


0.64 


0.75 


0.87 




+1579G>A 


0.61 


0.12 


0.35 


0.15 


0.93 


0.59 


0.47 


0.41 



Pvalues calculated from chi-square tests are presented for asthmatics with and without eosinophilia. Fvalues in boldface are significant (P<0.05). 
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Table 3. Gene-gene interactions between the CCR3 SNPs and eotaxin gene 
SNPs in blood eosinophilia in asthmatics 



Table 4. Analysis of gene-gene interactions by GMDR 



Model 


03 


SE 


Chi-square 


Pvalue 


totaxin- 1 oNr x LLnj etrect 










1010P^ A sy ppdi 
-\Jlob> AX LLnJeff 


n 1 1 
U. I I 


n /ic 
U.4b 


n nc 
U.Ub 


n Q1 
U.o I 


i7o a s/^v ppri 

-JZdA /UA LLnJeff 


-U.J I 


n 1A 
U.Z4 


1 70 

./u 


n 1 q 

U. I D 


i 1 11P^ A \y PPDO 

+ 1 ZJb> AX LLnJeff 


n co 
-U.bZ 


n oc 
U.Jb 


O 07 

J.U/ 


n no 
U.Uo 


i 171 A PPD1 

+3/ 1 l>AX LLnJeff 


-U. Id 


n oc 
U.Zb 


n oyi 
U.J4 


n cc 
U.Db 


+loo4b> AX LLnJeff 


n oc 
U.Zb 


n ai 
U.4Z 


0 /in 
U.4U 


n ro 
U.uo 


i 1C~71T^ PPD1 

+10/ 1 l>LX LLnJeff 


n nc 
U.Ub 


n oo 


n no 
U.UJ 


n 07 
U.O/ 


botaxin-Z bNrXLLhj ettect 










; ogp^* a \/ ppdi 
+ybb>AX LLnJeff 


U.U4 


n /io 
U.4Z 


n m 
U.U I 


n qo 

u.yj 


+304C>A(29L> l)X CCR3esf 


-0.68 


0.26 


6.89 


0.0087* 


+447C>TxCCR3 e u 


0.19 


0.23 


0.68 


0.41 


+1272A>GxCCR3efi 


-0.17 


0.35 


0.23 


0.63 


+1923A>CxCCR3 eff 


-0.14 


0.33 


0.17 


0.68 


Eotaxin-3SNPxCCR3 effect 










+77C>TxCCR3etf 


-0.10 


0.25 


0.15 


0.70 


+7164A>GxCCR3 ef f 


0.47 


0.34 


1.94 


0.16 


+1579G>AxCCR3 S ff 


-0.36 


0.46 


0.60 


0.44 



Model 



Accuracy Sign test (P) CVC 



Pvalues were calculated by logistic regression analysis. *Significant after the 
Bonferroni correction (P<0.01 for eotaxin-2). SE, standard error. 

174C> T and E0T2+304OA was the best model (accuracy= 
0.536, CVC 9/10; Table 4). This combination was also the best 
model after adjustment for age and sex as covariates. The com- 
bination of the CC genotype of CCR3-1 74C> T and AA or AC 
genotype of EOT2+304C>A was detected in the high-risk group. 
The combination of CT of CCR3-1 74C> 7/ and CC of EOT2+304C 
>A was also detected in the high-risk group (Figure). 

DISCUSSION 

Given the wealth of data showing that CCR3 and eotaxin par- 
ticipate in eosinophilic infiltration, the cooperative effect of the 
SNPs of CCR3 and eotaxin on eosinophilia in asthmatics was 
analyzed in this study. Three SNPs {CCR3-22557G>A, CCR3- 
520T>G, and CCR3-1 74C> T) in CCR3 were significantly asso- 
ciated with the number of eosinophils in patients with asthma. 11 
Associations between EOT2+1272A> G and the risk of asthma, 
and between EOTl+123G>A and high serum total IgE levels, 
were reported in Korean asthmatics. 15 However, these SNPs were 
not found as a best single-locus model in relation to blood eo- 
sinophilia in this study. 

Gene-gene interactions in various diseases have been studied 
using a variety of approaches. These studies demonstrated that 
the overall disease risk can be modeled as the product of the 
risk conferred by many independent risk factors. 16 In this study, 
gene-gene interactions were examined using 3 methods. First, 
a Chi-square test was conducted. Using this approach, EOT2+ 
304C>A was found to be significantly associated with the 3 in- 



CCR3-U40 T 

CCR3-1740 1 E0T2+304OA 
CCR3-520T> G, E0T1+361T>A, 

E0T2+304OA 
E0T1+3B1T>A, E0T2+304C>A, 

E0T2+447O I E0T2+1272A> G 



0.505(0.476) 0.62(0.828) 7/10 

0.536(0.549) 0.005(0.011) 9/10 

0.559(0.544) 0.011(0.011) 4/10 

0.525(0.576) 0.172(0.001) 3/10 



The numbers in parentheses are the results after adjustment for age and sex as 
covariates. Age was used as a continuous variable and sex as a discrete vari- 
able. CVC, cross-validation consistency. 

CCR3-174 OT 
C CT T 



A 




Accuracy 0.5355 
CVC 10/10 
P= 0.005 
ID'I=0.1871 



Figure. The best model was composed of CCR3-174C> Tand E0T2+304C> A. 
In each cell, the left bar represents a positive score and the right bar a negative 
score. High-risk genotype combinations are shaded dark grey, while low-risk 
genotypes are shaded light grey. Each P-value was obtained from a sign test, 
CVC, cross-validation consistency; |D'|: Lewontin's |D'| for linkage disequilibrium. 

tronic SNPs of CCR3 in patients with eosinophilia (Table 2). Next, 
we used logistic regression analysis and designed an epistatic 
model as described previously. 22 Three of the 4 CCR3 SNPs con- 
tributed independendy to the eosinophilic phenotype in asth- 
ma. This model can be used to evaluate the mean effect of all 
SNPs in CCR3. Of the 14 SNPs in the eotaxin gene family, only 1 
in eotaxin-2 {E0T2+304OA [29L>I]) was significandy associ- 
ated with blood eosinophilia; however, the association was 
modified by the effect of CCR3 SNPs {P= 0.0087; Table 3). Final- 
ly, we constructed an epistasis model using GMDR with and 
without covariate adjustment. Because the number of subjects 
in the 2 groups (asthmatics with or without eosinophilia) was 
different, GMDR was used instead of regular MDR. The combi- 
nation of EOT2+304C>A and CCR3-1 74C> 7/was the best mod- 
el. This model was still the best after adjustment for age and sex 
as covariates (Table 4). We successfully identified genotype 
combinations that contribute to eosinophilia in asthma using 
this model. Asthmatics, who had the C genotype of CCR3- 
1 74C> T and A or AC of EOT2+3040A and who had the CT 
genotype of CCR3-1 74C> T and C of E0T2+304OA (Figure), 
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were included in the high-risk group. Thus, the number of 
blood eosinophils is higher in asthmatics bearing these geno- 
type combinations. 

No single SNP in the eotaxin gene family was associated with 
blood eosinophilia in asthma (data not shown). However, 
EOT2+304C>A in which leucine is changed to isoleucine inter- 
acted significandy with the effect of CCR3 SNPs on blood eo- 
sinophilia in asthmatics. This finding remained significant after 
correction with the Bonferroni test (P<0.01). The eotaxin and 
CCR3 pathways are thought to interact biologically, suggesting 
a gene-gene interaction, or epistasis. The term "epistasis" has 
been used to describe a "non-independence of effect'.' In essence, 
if the effect of one unit is not predictable unless the value of an- 
other unit is known, the effect is epistatic. 25 In this regard, the 
effect of E0T2+304OA may not be disclosed if the effect of 
CCR3 SNPs is not considered. Moreover, there is no definite ev- 
idence that gene-gene interactions can be identified only among 
loci that show a significant association with phenotype. 

There are some limitations to this study. First, blood eosino- 
philia was included in our models instead of sputum eosino- 
philia; however, a small number of patients had sputum eosin- 
ophilia. Second, it is difficult to confirm the results of interac- 
tion models biologically. Eosinophils are further decreased in 
the lungs and peribronchial tissues of eotaxin-l/-2 double- 
knockout mice, compared with CCR3-deficient mice. 26 Howev- 
er, the knockout of eotaxin and CCR3 in animals has not been 
reported. 

In summary, novel gene-gene interactions between SNPs of 
CCR3 and an SNP of eotaxin-2 in the same phenotype were 
identified in this study. Our analysis of the epistatic influence of 
CCR3 and eotaxin gene variants suggests that these variants 
may be candidate markers for eosinophilia in asthma and could 
be important in understanding the genetics underlying allergic 
diseases. 
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